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Magnetoresistance (MR) has attracted tremendous attention for possible technological applications. 
Understanding the role of magnetism in manipulating MR may in turn steer the searching for new 
applicable MR materials. Here we show that antiferromagnetic (AFM) GdSi metal displays an anisotropic 
positive MR value (PMRV), up to —415%, accompanied by a large negative thermal volume expansion 

(NTVE). Around T N the PMRV translates to negative, down to 10.5%. Their theory-breaking 

magnetic-field dependencies [PMRV: dominantly linear; negative MR value (NMRV): quadratic] and the 
unusual NTVE indicate that PMRV is induced by the formation of magnetic polarons in 5 d bands, whereas 
NMRV is possibly due to abated electron-spin scattering resulting from magnetic-field-aligned local 4f 
spins. Our results may open up a new avenue of searching for giant MR materials by suppressing the AFM 
transition temperature, opposite the case in manganites, and provide a promising approach to novel 
magnetic and electric devices. 



Magnetoresistance (MR), a change in electrical resistivity when an external magnetic field (^ 0 ^) applied, 
occurs in metals, inorganic and organic semiconductors, and particularly close to an intermediate 
regime of the transformations of charge (insulator, metal) order and spin [paramagnetic (PM), fer- 
romagnetic (FM)] order in thin manganites' films as a colossal negative MR value (NMRV) 1 " 5 . The colossal MR 
(CMR) in manganites is accompanied by a shift of the transition temperature to a higher value by applied 
magnetic field so that a sharp peak appears in the MR values near the transition, e.g., the NMRV reaches 

95% at fi 0 H — 15 T near T c = 240 K in La 0 .85Sr 0 .i5MnO 3 6 . Although the double exchange interaction can 

qualitatively explain the CMR effect based only on the charge and spin degrees of freedom, the exact picture still 
remains elusive. 

We investigate here the intermetallic GdSi that crystallizes in the FeB-type structure (Pnma) and orders 
antiferromagnetically below —55 K 7 " 11 . The vanished orbital momentum (L = 0) naturally removes the crystal 
electric field (CEF) as well as its perturbation on magnetic interactions, which renders GdSi ideal for studying 
pure spin-magnetism-tuned transport. Meanwhile, there keep solely two main sources for the magnetovolume 
(MV) effect. One contribution is from the volume dependence of the Ruderman-Kittel-Kasuya-Yosida (RKKY) 
exchanges, which normally produces a negative MV (NMV) effect. The other is due to the spin-polarized 
itinerant moments in the conduction bands, resulting in a positive MV (PMV) effect. Therefore, the volume 
variation (NMV or PMV) with temperature in a magnetic state could be useful in determining the nature of its 
magnetic origin. Even after nearly 50 years of research on GdSi 8 " 11 , to our knowledge, there is little data on its 
temperature-dependent structural modifications and MR property. Here we report on detailed temperature- 
dependent powder diffraction, angular- dependent magnetic characteristic and transport studies of GdSi single 
crystals. We find anomalous anisotropic giant MR and spontaneous magnetostriction (MS) effects, in particular, 
an antiferromagnetic (AFM) -driven negative thermal volume expansion (NTVE) and a nontrivial positive to 
negative MR transition, which can be well understood by combining the magnetic tunnel of conduction electrons 
and the concept of magnetic polarons 12 , i.e., local short-range FM spin regimes 12 . 
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Figure 1 | Temperature-dependent structural parameters. 

(a) Anisotropic character of the lattice-constants, a, b and c, variation. 

(b) Anomalous unit-cell volume, V, expansion with temperature in the 
Pnma symmetry. The solid lines are theoretical estimates of the 
temperature-dependent structural parameters using the Griineisen model 
with Debye temperature of 0 D = 340 K, which is the same as reported in 
Ref. [10]. Upon warming, two appreciable anomalies display in the 
structural parameters at respective temperatures of T N1 ~ 54.3 K (at 
0.06 T) and T\ ~ 100 K. Error bars in (a) and (b) are the standard 
deviation obtained from the Fullprof refinements. 

Results 

Structural studies. The x-ray powder diffraction analysis (Fig. 1) 
clearly shows three distinct structural regimes (I, II and III). Upon 
cooling, the refined (Re) a, b, cand V shrink almost linearly before T l5 
followed by a slower decrease until sharp turns at T N1 , an onset 
temperature of the AFM transition. Below T N1 , the increases of a 
and b and the decrease of c ultimately result in an unusual NTVE (i.e., 
PMV) in the unit- cell volume V, which is quite useful in fabricating 
applicable materials with controlled thermal expansion values 1314 . 
Obviously broadening in the nuclear Bragg peaks is present in the 
I regime (see Supplementary Fig. Sib), which is attributed to the 
magnetoelastic effect (the coupling between magnetic moments 
and lattice strains). The strain distribution patterns are extracted 
and shown in Supplementary Fig. S2. The variations in a, b, c and 
V below T N1 imply that magnetic anisotropy (MA) (which is 
consistent with Supplementary Fig. S3a and uncommon for the S- 
state Gd- compounds), spontaneous PMV and anisotropic MS effects 
exist in GdSi. 

The CEF is mainly responsible for the giant MS effect in rare-earth 
(RE) compounds. This effect is thus expected to be negligible in 
GdSi. However, below 7\, structural parameters shown in Fig. la 
obviously deviate from the theoretical estimates (solid lines) by the 

a 10K — a 10K foioK_uoK 
Gnineisen (Gr) law, e.g., Re 1QK Gr =4.94(6) x 10~ 4 , Re Gr 

a Gr ^Gr 
C 10K _ C 10K 

= 8.24(9) x 10- 4 and Re 1QK Gr =-2.60(8) x 10 " 4 , denoting large 

C Gr 

anisotropic spontaneous MS effects. The formation of long-range- 
ordered (LRO) AFM state is a process of the growth of sublattice 



FM domains. The enlargement of FM domain volumes with decreas- 
ing temperature may accumulate strains on the domain walls, which 
is the microscopic mechanism for the magnetic-field-induced MS 
effect in ferromagnets. Therefore, including the effect of the molecu- 
lar field of one Gd- AFM- sublattice on the other is indispensable to 
understand the spontaneous MS effect in GdSi. According to the 
Stoner model for itinerant magnetic electrons, the positive magnetic 
pressure P M associated with the magnetic ordering in a band is 

proportional to ^^^ M 2 , where D, V and M represent the elec- 
tronic density of states at the Fermi energy, the volume and the 
magnetic moment, respectively. The spontaneous PMV effect (i.e., 

ylOK _ ylOK 

NTVE), e.g., Re 1QK Gr = 1.04(2) x 10" 3 , is therefore attributed 

^Gr 

mainly to the increases of D (corresponding to the pronounced 
decrease of p below T N in Fig. 2a) and the induced itinerate spin- 
moments in conduction bands. Similar MS and MV effects were also 
reported in Gd 3 Ni 15 , where, however, they are ascribed to the itin- 
erant character of the Ni 3d bands. 
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Figure 2 | Temperature variations of resistivity and MR effect. 

(a) Resistivity measurements with current I along the a and c axes under 
applied magnetic fields of 0, 1 and 8 T. T\ = 100 K labels the temperature 
where one structural anomaly occurs as shown in Fig. 1 . The dashed lines are 
fits between 10 K and 40 K (details in text), and extrapolated to higher 
temperatures, (b) Corresponding MR values versus temperature. The MR 
effect along the a axis has a similar trend to that of the c axis albeit with a 
lower value. The positive MR values decrease sharply with increasing 
temperature below —20 K, then gradually transfer into negative around T N , 
and persist up to ~ 120 K. By contrast, in amorphous Gd^Si^ films, only the 
NMRV was observed 16 " 18 near the metal-insulator transition. In (a), there is 
no big difference for the data at 0 and 1 T above —20 K. The solid lines in (a) 
and (b) in the dashed arrow direction in turn correspond to the ordinal axis- 
direction and applied magnetic-field strength as labeled, respectively. 
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Resistivity measurements. As depicted in Fig. 2a, the zero -field- 
cooling (ZFC) electrical resistivity at 0 T with current I along the a 
and c axes decreases linearly due to the weakened thermal excitations 
upon cooling until around T N (vertical bar), where a hump with 
negative slop appears along only the c axis, probably attributed to 
the magnetic superzone effect 19 as a consequence of the AFM order- 
ing. Below T N , they decrease steeply like a reasonable AFM metal. 
The resistivity between 10 K and 40 K 20 can be well fit to p(T) = p 0 + 
kT™, shown as dashed lines. This produces w c0T = 1.51(2), w c,8T = 
1.56(4), w a ' 0T = 1.43(2) and w fl ' 8T = 1.51(1). All w values are much 
smaller than 5 20 indicative of anisotropic magnetic interactions. 

Anisotropic MR effect and positive MR value (PMRV) to NMRV 
transition. The most intriguing results from resistivity measure- 
ments are the anisotropic MR effect (Figs 2b, 3 and 4) and the 
existence of both positive and negative MR values, up to —415% 
(comparable to the CMR value in manganites 1,6 and one to two 
orders of magnitude larger than that of the RE-metals 21 ) and down 

to 10.5% along the c axis at 8 T and 3 K and 52.8 K, respectively. 

The MR anisotropy in the ac and be plans is shown in Fig. 3. They 
display a twofold symmetry at 7K (Figs 3a and 3d). We notice that 
applied magnetic field of 8 T does not suppress (produce) the (a) 
hump along the c and a axes, respectively, near T N , and the MR 
twofold symmetry is persistent from 1 to 8 T, indicating that app- 
lied magnetic field in a strength of 8 T may just align or localize the 
5d moments, and slightly rotate and tilt the 4f moments while 
conserving the superzone energy gap. Therefore, the MR effect in 
GdSi exhibits a well separate feature of the temperature regions for 



the positive and the negative MR values (Fig. 2b), respectively, which 
is induced jointly by the AFM superzone effect 19 and the shift of the 
AFM transition to lower temperatures in external applied magnetic 
field analogous to the case in manganites 1,2,6 . For metals, mean-field 
theories predict that spin fluctuations induced by applied magne- 
tic field from the localized magnetism produce a PMRV with the 
quadratic-field dependence in antiferromagnets, whereas a NMRV 
with the linear variation in ferromagnets and paramagnets 22 . 
However, in GdSi, the PMRV in the AFM state mainly displays a 
linear magnetic-field dependence (Figs 4a and 4c), while above T N1 
the absolute NMRV is proportional to the square of the strength of 
applied magnetic field (Figs 4b and 4d). Both the positive and 
negative MR effects do not saturate at utilized maximum fi 0 H = 
8 T (Fig. 4). In addition, the ratio of the resistivity at 160 K and 
7 K in Fig. 2a is already —12-25, therefore, the cyclotron motion 
of the conduction electrons could be neglected at /n 0 H = 8 T (i.e., 
co c t<C1, where w c is the cyclotron frequency and t is the life time of 
the conduction electrons). Therefore, these uncommon magnetic- 
field variations indicate new transport mechanisms for the MR 
effects of GdSi. 

Discussion 

In GdSi, the conduction electrons (mainly 5d bands) are different 
from those responsible for the magnetism (4/ component plus pos- 
sible part of the 5d component). The former is normally delocalized, 
acting as the magnetic glue among magnetic ions (Fig. 5), and scat- 
tered by them, leading to electrical resistance. The magnetism from 
the 4/ part is generally localized with weak interactions. Therefore, 
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Figure 3 | Angular-dependent MR values under different applied magnetic fields, (a) Angular-dependent MR values with applied magnetic field in a 
range of 0 to 8 T (1 T step) at 7 K, and (b) at 120 K (0, 2 and 8 T). (c) For (a) and (b) measurements, current lis along the crystallographic c axis, and 
applied magnetic field, /UoH y rotates away from the a axis with an angle of 0 in the ac plane, (d) Angular-dependent MR values with applied magnetic field 
in a range of 0 to 8 T(l Tstep) at 7 K, and (e) at 120 K (0,2 and 8 T). (f) For (d) and (e) measurements, III a-axis and applied magnetic field, fi 0 H, rotates 
in the be plane with an angle of y deviated from the b axis. 
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Figure 4 | Field- and temperature- dependent MR values, (a) Field and temperature dependencies of the MR values with I\\ c-axis and applied magnetic 
field, fi 0 H y along the a or c axis below T N , and (b) above T N . The representatives of the linear-field dependence of the PMRV (i.e., the PMRV is 
proportional to the strength of applied magnetic field) at 7 K (blow T N ) and the quadratic variation of the NMRV (i.e., the absolute NMRV is 
proportional to the square of the strength of applied magnetic field) at 60 K (above T N ) were shown in (a) and (b), respectively, (c) Field and temperature 
dependencies of the MR values with I\\ a- axis and applied magnetic field, /u 0 H y along the b or c axis below T N , and (d) above T N . 



the LRO AFM state originates mainly from the isotropic RKKY 
interactions through conduction bands 23 . The interaction between 
localized moments, M\ oc , and itinerant ones, m iti , can generate an 
extraordinarily large Zeeman splitting in the mean-field approxi- 
mation 24 



E=g*fi B m i ti H + 2/ (H) m i ti ■ ( M L 



(1) 



where g* is the spectroscopic splitting factors for the carriers, fi B is 
the Bohr magneton, / is the effective exchange coefficient and 
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Figure 5 | Schematic illustration of the spin states with and without applied magnetic field, /£ 0 H, in different temperature regimes, (a-c) At zero 
magnetic field, (d-f) At applied magnetic field of /u 0 H. When T > T l5 spin moments more or less rotate (f), depending on the strength of /u 0 H and the size 
of MA, from a pure PA state (that is strictly observing the Curie-Weiss law as shown in Supplementary Fig. S4) in (c). When T N1 < T< T\, the short-range 
AFM spins that are attributed only to the local 4/moments appear (b) accompanied by the generations of polarized itinerate 5d spins, based on the 
deviation of the unit-cell volume from the Griineisen model shown in Fig. lb, and possible small amount of localized 5d spins according to equation ( 1 ). 
Applied magnetic field mainly aligns the local AFM spins (e), leading to a decrease of the electron-spin scattering and resultant the NMRV. When T < 
T N2 , the LRO AFM state with almost equivalent AFM and FM interactions (see Supplementary Fig. S4) forms (a) with more itinerate 5d moments (based 
on the large PMV effect shown in Fig. lb). Applied magnetic field mainly localizes more 5d moments by enhancing the exchange of /in equation (1), 
resulting in the formation of magnetic polarons and the consequent PMRV (d). 
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(Mioc) is the averaged local moment in the regime of band electrons. 
Since m iti is usually small, and the second term could be very large 
(e.g., amounting to fractions of an eV in the LRO magnetic state of 
Eu- compounds 24 ), in addition, / is strongly associated with applied 

magnetic field by virtue of modifying spin fluctuations of M\ OC) the 
formation of magnetic polarons in the 5d bands by this splitting in 
the LRO AFM state of GdSi is thus possible. Therefore, when T < 
T N2 , the modified /at 8 T drives some of the conducting moments (as 
foregoing remarks) to form local magnetic polarons that lead to a 
largely degenerate conduction (i.e., PMRV) (Figs 5a and 5b). In this 
case, the more extended 5d bands almost certainly offer a small FM 
component, which dominates the linear-magnetic-field dependence 
below T N2 . 

Above T N1 , the LRO AFM 4/moments disappear, which is implied 
by the change of the x relationship between the a, b and c axes (see 
Supplementary Fig. S3a). However, the upward deviation of 1/x from 
the Curie- Weiss law below T x (see Supplementary Fig. S4) along with 
the concomitant anomalous volume lattice distortions (Fig. 1) indi- 
cates that a kind of striking AFM state exists there, which can be 
ascribed only to the local short-range AFM 4f moments (Fig. 5b) 
according to equation (1) (broad 5d states indeed could not form a 
magnetic state by themselves). This local AFM state theoretically 22 
leads to the quadratic-magnetic-field variation (Figs 4b and 4d). The 
local AFM moments above T N1 are possibly oriented irregularly in 
the PM state and become compulsorily aligned and to some extent 
spin-tilted (Fig. 5e), leading to the decrease in resistivity and thereby 
promoting a NMRV. 

Including the contribution of magnetic polarons into the mean- 
field approximation of the MR effect of AFM metals 22,25 " 27 , the MR 
value (MRV) then equals 

|MRV|~A|i< f + M^| 2 

( 2 ) 

/i^4f| 2 i—ylfll— >5d| I >5d | 2 \ V ' 

~A(jM In | +2|M In ||M In |cosa+|M In | J, 
where A is a constant, M In is the induced local magnetization by 

_»4f -^5d 

applied magnetic field, /n 0 H, a is the angle between M In and M In , 

M In = x inMo^' an d X in * s me magnetic susceptibility of the 4/ site 
along the AFM sublattice direction. At T N1 < T < T l5 neglecting 

l^ 5d l l^ 4f | 2 I ^4f 1 2 

M In , | NMRV | gc M In = \Xin\ rt)-" > resulting in a quadratic- 
magnetic-field dependence. At T < T N2 , in our case, the linear part 

I _^4f I I -^5d I I _^ 4f I 

in equation (2) is dominate, i.e., PMRVoc M In M In cos a= X In 

I — >5d I 

M In Lu 0 Hcos a, which produces not only the linear-magnetic-field 

dependence but also the twofold MR symmetry (Figs 3a and 3d). 
Quantitative analysis requires the knowledge of the exact directions 

_^ 4f _^5d 

of x in an d Afin • The involvement of magnetic polarons in under- 
standing the MR effects here is further supported by the existence of 
strong FM Gd-Gd interactions (see Supplementary Fig. S4). 

RE 4f electrons generally remain localized so that their properties 
in an alloy closely resemble those in the single-elemental metals. It is 
therefore interesting to compare the microscopic conducting mech- 
anism of GdSi with that of the 3d transition metals 28 ' 29 . Roughly, the 
4f shells highly screened by the 5s and 5p states are better shielded 
than the transition metal (TM) 3d ones. Therefore, the 4f electrons 
are well embedded within the atom, and the 5d and 6s states act as 
conduction electrons in the metals, though band structure calcula- 
tions suggested that the upper part of the 4f spin minority may 
hybridize with the unoccupied 5d and 6s bands 30 . On the contrary, 
for 3d TMs, the conductivity is mainly dominated by the s-s and d-d 
transitions at low temperatures, in addition to the s-d transitions 
(increasing the effective electron mass and shortening the mean free 



path) at high temperatures 28,29 , supposing that both s and d electrons 
are conduction electrons. Since almost all RE 4/metals and only some 
of the 3d TMs (Fe, Co, Ni and Mn) order magnetically at low tem- 
peratures, magnetism plays a more crucial role in controlling elec- 
trical property of RE 4f metals. 

In summary, GdSi displays a wide array of novel behaviors includ- 
ing a two-step AFM transition, an anisotropic spontaneous MS 
effect, magnetostrictive strains, in particular, an exotic NTVE (spon- 
taneous PMV) and a peculiar positive to negative MR transition, 
indicating a strong coupling between spin, lattice and charge degrees 
of freedom. The peculiar PMRV to NMRV transition and their non- 
trivial 22 magnetic-field dependencies are predominated by the nature 
of the 4/moments and the polaronic 5d carriers. Coexistence of these 
extraordinary behaviors in the same Gd-compound is unique. The 
present results make GdSi an fascinating system for theoretical and 
further experimental explorations. 

Methods 

Since the neutral Gd is a strong thermal neutron absorber, we therefore performed a 
powder x-ray diffraction study of the pulverized GdSi single crystal on an in-house 
diffractometer employing the copper K al = 1.5406(9) A as the radiation with a 29 
step size of 0.005° from 10 to 300 K to explore the temperature-dependent structural 
modifications. X-ray powder diffraction data were analyzed by Fullprof suite 31 . ZFC 
dc magnetization (M) measurements were performed on a Quantum Design MPMS- 
7 superconducting quantum interference device (SQUID) magnetometer. The ZFC 
electrical resistivity (p) of bar-shaped (~ 0.6 X 1 X 6 mm 3 ) single crystals was 
measured by standard DC four-probe technique using a commercial physical prop- 
erty measurement system (PPMS), equipped with the option of rotating the sample in 
situ. Commercial silver paint and 25 fim gold wire were used for electrical contacts. 

p{0,T,(p/y) 

p(0, T, <j)/y) are the resistivity with and without fi 0 H at a given temperature, T, and a 
rotating angle of (f) or y (Figs 3c and 3f). 
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